Geometric transition in friction for flow over a bubble mattress 
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Laminar flow over a bubble mattress is expected to experience a significant reduction in friction since the 
individual surfaces of the bubbles are shear-free. However, if the bubbles are sufficiently curved, their protrusion 
into the fluid and along the flow direction can lead to an increase in friction as was recently demonstrated 
experimentally and computationally. We provide in this paper a simple model for this result. We consider 
a shear flow at low Reynolds number past a two-dimensional array of bubbles, and calculate analytically the 
effective slip length of the surface as function of the bubble geometry in the dilute limit. Our model is able to 
reproduce quantitatively the relationship between effective friction and bubble geometry obtained in numerical 
computations, and in particular: (a) The asymmetry in friction between convex and concave bubbles, and (b) 
the existence of a geometric transition from reduced to enhanced friction at a critical bubble protrusion angle. 



Microfluidic devices are used to manipulate small volumes 
of liquids (micro- and nano-liters), and have important appli- 
cations in biology, chemistry and engineering |[lll2l • However, 
from a practical standpoint, one important caveat of decreas- 
ing length scales is the increase in friction. For flow driven by 
a pressure gradient, Vf , in a device of typical cross-sectional 
scale a, the volume flow rate scales as Q ~ a'VPIrj, where 
rj is the viscosity of the fluid. Consequently, imposing a con- 
stant flow rate in a device decreasing in size requires a sharp 
increase of the applied pressure gradient as VP ~ l/a^. For 
example, driving a flow rate of one microliter of water per sec- 
ond in a 10 jum-wide and 1 cm-long device requires a pressure 
drop of many atmospheres. 

Viscous friction is therefore an important issue on small 
scales, and engineers have looked for ways to reduce it. Re- 
cently a number of groups have reported slip, or apparent slip, 
for simple shear or pressure-driven flows over hydrophobic 
surfaces (see Refs. ij, 0, Htl for a review). Slip is typically 
characterized by a slip length. A, which is the fictional dis- 
tance below the surface at which the velocity field would ex- 
trapolate to zero. For a shear flow with shear rate y over a flat 
slipping surface, the velocity field is written as 

u^(Ay + fy)x (1) 

where x is the flow direction, and y the direction perpendic- 
ular to the surface. For a system of typical cross-sectional 
size a, friction decreases as some increasing function of A/a 
fi 0j S la, S S IM] • Therefore for small enough systems, or 
large enough slip lengths, the friction reduction could be sub- 
stantial. However, intrinsic slip lengths of hydrophobic sur- 
faces do not exceed tens of nanometers ||3l|j,|5|l- As a result, 
the effect is small for devices on the micron-scale, and other 
means of friction-reduction must be devised. 

One idea recently discussed is to exploit super-hydrophobic 
surfaces 1 10, 1 1] to reduce drag. When an hydrophobic solid 
surface (Fig.[T^) is sufficiently rough, under certain geometric 
conditions, a low-pressure high-surface energy fluid in contact 
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with the solid can spontaneously de-wet [11], thus transition- 
ing from a Wenzel state where the fluid fills the grooves on the 
surface (Fig. \^), to a fakir-like Cassie state where the fluid 
sits part ially on the solid surface and partially on gas (Fig.[T]:) 
ifToi lull. This phenomenon can also occur at the nanometer 
scale uM- When de-wetting does occur, shear stresses on the 
portion of the fluid interface in contact with the gas are ex- 
pected to be insignificant, and the effective friction of such a 
super-hydrophobic surface should be significantly lower than 
that of a solid surface. Drag reduction by super-hydrophobic 
surfaces was indeed demonstrated experimentally, both for 
macro-scale pipe flow 1 13], as well as microchannels fT^IHtl- 
Numerical simulations also show a similar effect at the molec- 
ular scale I l3l ■ 

Liquids in contact with a mixed solid/gas interface are 
therefore expected to display low friction. However, one issue 
can potentially limit the drag-reducing properties of such sur- 
faces, namely the geometry of the gas-liquid interface. When 
the pressure in the gas is similar to that in the liquid above it, 
the interface remains flat (Fig.[T]l), and drag is expected to be 
reduced by the gas-Uquid interface. However, when the pres- 
sure in the gas is larger than that in the liquid, the interface 
becomes curved (Fig. [1^), with mean curvature given by the 
Young-Laplace equation lITlll . In that case, the protrusion of 
the interface into the liquid curtails the drag-reducing prop- 
erties of the surface by distorting the flow streamlines in the 
shear direction HI O [H, [H] . Notably, if we denote by 6 




FIG. 1: Surperhydrophobic surfaces and bubble mattress: Liquid 
drop on a: (a) flat hydrophobic surface; (b) rough hydrophobic sur- 
face in Wenzel state; (c) rough hydrophobic surface in Cassie state. 
When the drops sits partially on air (Cassie state), the gas-liquid in- 
terface can be flat (d) or curved (e), with a protrusion angle 6 mea- 
sured with respect to the horizontal direction. 
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FIG. 2: Setup for the two-dimensional calculation and notation: (a) 
shear flow (shear rate, y) over a collection of bubbles of projected 
radius c, protrusion angle 9 and separation distance 2L. The calcu- 
lation is performed asymptotically in the dilute limit, (f> = c/L <k 1; 
(b): The case 6 > corresponds to a concave bubble protruding into 
the liquid; (c): 6 < models a convex bubble curved away from the 
otherwise flat surface. 



the value of the protrusion angle of the spherical-cap interface 
with respect to the horizontal surface (see Fig.[T^), two recent 
experimental and numerical studies demonstrated that if the 
interface extends into the liquid beyond a critical value Oc, the 
friction of the surface becomes larger than that of the origi- 
nal (flat) no-slip surface | 8„,S1- These results imply therefore 
that the geometry of free surfaces plays a crucial role vis-a-vis 
the resistance to viscous flow, implying in particular that the 
recently-discovered nano-bubbles 1.20. 21] might not always 
lead to as low a friction as previously thought. 

Here, we specifically address this geometric transition from 
reduced to enhanced friction for shear flow over a collection of 
bubbles. We consider a simple two-dimensional model to pre- 
dict the critical value of the apparent angle 6c of the bubbles 
on the (otherwise) horizontal surface at which the effective 
surface starts displaying a friction larger than that obtained by 
a smooth solid surface. In the two-dimensional case, the ve- 
locity field can be solved exactly in the dilute limit, and the 
geometric transition is seen to happen at Oc ~ 65°. More gen- 
erally the relationship between effective slip length and ap- 
parent angle from numerical simulations is well represented 
by our simple model |1^ |^ . 

We consider therefore the setup displayed in Fig. |2l A di- 
lute collection of two-dimensional bubbles (no-shear surfaces) 
on an otherwise no-slip surface. We assume the capillary 
number to be sufficiently small that the bubbles have perfectly 
circular shape. Their projected radius on the surface is c and 
they are separated by a distance 2L. The bubble surface cov- 
erage is therefore given hy (p - c/L. The protrusion angle of 
the the bubble into the fluid is denoted 0; > is a concave 
bubble that protrudes into the fluid (Fig.|2]3), whereas < is 
a convex bubble protruding underneath the surface (Fig. |2j;). 
The velocity field away from the surface is a pure shear flow, 
with constant shear rate y. We denote by x the direction along 



the surface, y the direction of shear, and z the third direction, 
and solve for the flow field in the limit of low Reynolds num- 
ber il. 

We first examine the case of a single bubble. The ve- 
locity field, V, is solved by introducing a streamfunction, 
V - curl(i^z), where W^ifr = in the half-plane y > 0, indented 
by an arc of the circle + (y + ccotO)^ - c^csc^O which per- 
turbs the shear flow i/r = \yy^- There is perfect slip without 
penetration at the surface of the two-dimensional bubble, and 
no slip is enforced at the rigid plane y - Q. We introduce 
toroidal coordinates defined by 

c sinh ^ 



cosh^ + cos 77' 



y = 



csin?7 
cosh^ + cos 77' 



(2) 



which is the conformal mapping x + iy - c tanh(^ -1- irf) / 2 . The 
fluid region is then -00 < ^ < 00,6 < t] < n, with ^ - +00 
at the special points (+c, 0), ^ = on the y-axis, rj - 6 on the 
arc and 77 = ;7r at 37 = 0, > c. A suitable biharmonic stream 
function, asymptotic at infinity to the prescribed shear flow, is 
given by 1231 12411 



ip — yc 



1 • 2 
5 sin 77 



fis,r]) cosi^di 



(cosh^ + cos 77)2 
sinh s(n - 77) 



-A(s) sin?7 



+B(s) 



cos 77- 
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sinh i(7r ■ 



1) 



cosh^ + cos 77 



+ sin 77 cosh s(7t - 77) 



(3) 



Here, two functions have been eliminated by imposing no slip 
(i/r = = di/z/dri atrj-n) on the rigid plane on either side of 
the bubble. Note that using the method in Eq. (O allows us to 
capture the appropriate stress singularity at the contact point 
between the bubble and the surface. 

The first equation for the functions A and B is established 
by demanding that the arc 77 = be a streamline on which 
tp -0. Thus the transform. 
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sinh S77 
sinh S7T 



cos ds - 



sin 77 



(cosh^ + cos 77) ' 



yields, by substituting Eq. (O, 

f(s, 6) = sine- 



sinh s6 



(\r]\<n), (4) 



(5) 



sinh S7T 

The second equation is established by demanding that the 
strain rate component, e^^i, vanishes at 77 - in order to 
achieve perfect slip. Since 

d I I dtff 



drj \ drj 



cosh^ + cos 77' 
(6) 



it is found that 
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(9) 
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The first and second derivatives of Eq. (|4|i give 

2s[s sin 0sinh s6 + cos ^cosh s6] 



sinh S7T 



(10) 

After substitution of Eq. (O into Eq. Q and Eq. ( fTOl i. it is 
found that 



A(s) 



sinh 2s(7T - 6) + s sin 26 

5 sin 20 cosh sn: + sinh s{n - 29) 



(11) 



and 



cos 26 ■ 



B(s) 



sinh sn 



s sin 26 



sinh 2s{n - 6) + s sin 26 



(12) 



We can now derive the resuh for a period array of bubbles. 
Away from the surface, the far-field velocity field is given, 
according to Eq. Q, by ^'^ + {n - 77)^ — » (// < tt), with 



2c^ 



y ' 



2c{n — rj) 



(13) 



Hence the asymptotic form of the perturbation term in Eq. (|3]l 
is given by 



2y& 



y 
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A{s)ds. 



(14) 



In the presence of a periodic array of bubbles, with period 2L, 

„2 



the prefactor in Eq. ( fT4l l is replaced by > 



y 



(x - 2nLy + 



with a mean value available from the identity 

y^dx y r°° ydx 



-\f- 



- 

CO ^ 



2nL)2 + 3;2 

Consequently, the far field velocity perturbation is given by 



nyc 



which, with A(s) given by Eq. (fTTI ). is identified as an effective 
slip speed at y = with dimensionless slip length (see Eq.[T]i. 



A Ic 

— -TlX — 

c \L 



A{s)ds. 



(17) 



Let us now compare the predictions of our model (Eq. [TTJ 
with the numerical computations of Refs I's', The results 
are displayed in Fig. [3] where we plot the effective slip length 
of the surface (non-dimensionalized by the projected bubble 
radius on the surface, c) as a function of the protrusion angle, 
6, into the fluid. The simulations of Ref. [8] consider a three- 
dimensional square lattice, with a surface coverage of the bub- 
bles (p - 0.68 (squares). The computations of Ref. |9] are 
for three different lattices (square, circles; rectangular lattice, 
lozenges; rhombic lattice, triangles) and in all cases where the 
dependence on the protrusion angle is studied the fraction of 
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FIG. 3: (color online) Comparison between the numerical data of 
Ref. (§] (empty symbols) and |9] (filled symbols) and our two- 
dimensional dilute model: Dimensionless effective slip length as a 
function of the protrusion angle of the bubble into the fluid. Sym- 
bols: Square lattice (□, <^ = 0.68 |8]; •, ^ = 0.43 Hi); rectangular 
lattice (♦, = 0.43 [9]); rhombic lattice (a, = 0.43 fP]). Lines: 
Our two-dimensional dilute model, Eq. for two different sur- 
face coverage, (f> = 0.43 ( ) and = 0.68 ( — ). 



the surface covered by the bubbles is = 0.43. We also plot 
in Fig. 13 the results of our model for both surface coverage 
considered in the two previous studies: cp = 0.43 (dashed line) 
and (f> = 0.68 (sohd fine) fM - 

The main features of the full numerical results are seen to 
be reproduced by our analytical model. There exists a criti- 
cal protrusion angle, 6^, above which the effect of the wall- 
attached bubbles displays a transition from reduced (6 < 6c) 
to enhanced friction (6 > ft ). Our model predicts ft ~ 65°, 
in good agreement with the results of Ref. |8] (ft- w 62°) and 
Ref. 1^ (ft ~ 69°). Our model also successfully reproduces 
the asymmetry in the friction curves between 6 <Q and 6 > 0, 
indicating a qualitative difference between the effect of con- 
vex and concave bubbles on shear flow. Experimentally, this 
asymmetry implies that flow over an array of bubbles at a 
lower pressure than the fluid (6 < 0) is less likely to show 
an increase of wall friction than flow over an array of bubbles 
with a larger pressure than the fluid {6 > 0). 

The model we consider here is two-dimensional. A simple 
argument can also be made to show that the critical protru- 
sion angle lies in the range 0° < ft. < 90° in three dimen- 
sions. For 6 = 0°, the bubbles are flat disks, for which the 
calculation is available in Ref. ll25ll . where it is found ll29tl that 
the asymptotic form of the perturbation flow component v ^ is 
2yc^yx^ /37Tr^ , with - x^ + y^ + z^- In the presence of a 
(dilute) square lattice of such disks, with period 2L, the ef- 
fective slip length is simply found to be /l/c - c^/9L^. This 
is the same result as that obtained by Ref. |@], and indicates 
that bubbles with 6-0° reduce the effective friction of the 
surface, and therefore 0° < ft. In the case where 6 = 90° 
(hemispherical bubble), the velocity field is readily verified to 
be given by v = - c-'jcf /r'^j ifTslfTgll . In the presence of a 
square lattice of such bubbles, the perturbation velocity in the 
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far field can be simply added up (in the dilute limit), and the 
obtained slip length is given by /l/c = -jic^jbL^. The negative 
slip length indicates that the array of stress-free hemispheres 
increases the effective friction of these surfaces, and therefore 
Be < 90°. 

In summary, this paper presents a two-dimensional model 
of shear flow past an array of bubbles. We calculate the ef- 
fective slip length of the surface in the dilute limit, and found 
good agreement with recent numerical three-dimensional nu- 
merical simulations for flow over super-hydrophobic surfaces. 
In particular, our results reproduce the asymmetry between 
convex and concave bubbles, as well as the existence of a ge- 
ometric transition from low to high friction at a critical bubble 
protrusion angle. Our approach provides therefore the min- 
imal model necessary to quantitatively capture the interplay 



between low friction and geometry in shear flows. 

Our results could be extended by relaxing the various as- 
sumptions made in the paper, in particular by going beyond 
the two-dimensional, dilute and zero Capillary number limits. 
In the case of three-dimensional bubbles, the issue of geome- 
try of the lattice needs to be further explored as well. Experi- 
mentally, the fabrication of controlled porous surfaces 112611 or 
pre-designed bubble lattices f2!f\ would present exciting op- 
portunities to map out and optimize the frictional properties 
of super-hydrophobic surfaces. In addition, simulations show- 
ing an increase in friction of a bubble mattress with capillary 
numbers |9] could be addressed with a similar framework. 
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